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ABSTRACT 



Method and apparatus for non-invasive testing of digital 
communications systems. Amplitude measurements are 
made for multiple frequencies of a multi-frequency com- 
munication system, converted to the time domain. An adap- 
tive filter output is matched to the time domain representa- 
tion to characterize the channel. Impedance mismatches may 
be precisely located using this technique. An error signal 
representing a difference between a signal transmitted 
through the channel and a received signal is estimated and 
analyzed. The error signal is separated into components 
corresponding to contributions by wide band noise, residual 
phase modulation, and residual amplitude modulation. Iden- 
tification and removal of narrow-band interferers may occur 
prior to this separation. Bit error rate and system margin 
computations employ a Monte Carlo simulation of the 
various error sources. This provides a well refined estimate 
of bit error rate and system margin. 

21 Claims, 16 Drawing Sheets 
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NON-INVASIVE DIGITAL CABLE TEST estimated and analyzed. The error signal is separated into 

SYSTEM components corresponding to contributions by wide band 

noise, residual phase modulation, and residual amplitude 

STATEMENT OF RELATED APPLICATIONS modulation. Identification and removal of narrow-band 

This application is a division of Ser. No. 08/791,386, filed 5 inte & K ' s ma y P 1 ™ to separation. Bit error rate 

Jan. 30, 1997, now VS. Pat. No. 6,061,393, which is a *? d ? y ^ tem 1 mare ? n comP"<ations employ a Monte Carlo 

continuation-in-part of U.S. patent application Ser. No. smu l an ° n sunulates the various error sources. This 

08/532,017, filed Sep. 21, 1995, now U.S. Pat. No. 5,757, P rovl ?es a well refined estimate of bit error rate and system 

766, which is in turn a continuation -in-part-of U.S. patent margin. 

application Ser. No. 08/430,909, filed Apr. 27, 1995, now 10 In accordance with a third aspect of the invention, a 

abandoned. Hie contents of both parent applications are cluster variance measurement based on an error signal is 

herein incorporated by reference for all purposes. used t0 estimate bit error rate. In accordance with a fourth 

aspect of the invention, an error signal is displayed in the 

BACKGROUND OF THE INVENTION frequency domain relative to a measured power of a 

_ 15 received signal. 

The present invention relates to measuring signal deera- ™ . t . M1 . , „ , 

j . * • . . . T . , , . to The invention will be better understood bv reference to 

dation in communications systems. In particular, the inven- tU ~ „ , , , , ... . v , , 

t; f _ . ma *u~A a ♦ f • - i tn e following detailed description m connection with the 

Uon relates to methods and apparatus for non-invasively • h ■ 

measuring performance of a digital broadcast system. accompanying awings. 

Digital communication technologies offer numerous 20 BRIEF DESCRIPTION OF THE DRAWINGS 

advantages over their analog predecessors. Nearly error-free FIG. 1 depicts a simplified representation of a represen- 

transmission is guaranteed if a threshold signal-to-noise tative digital communication system suitable for testing in 

ratio is available. Channel distortions may be corrected accordance with the present invention 

using adaptive equalizers Coding techniques may be advan- nG 2Ais a block ^ d ictin OQe im p lemeDtatioa 

tageously employed both to overcome channel-related sig- 25 of a ^ of a ril - - fal reccivf £ 

nal impairments and to minimize usage of bandwidth. . t „ , . , . „ . 

° FIG. 2B which depicts a so-called signal constellation for 

Nonetheless, digital commumcaaon system operation is a representative digital modulation scheme, 64 point QAM. 

not fault-free. The transmission of information over a physi- CT r ^ ^« w _ *■ * i * «■ c 

i j* . * . . / FIG. 2C depicts an alternative implementation of a por- 

cal medium ultimately requires analog components such as ^ e , . r ^ 

nr t-£ •« i m tion of a digital receiver. 

RF mixers, amplifiers, oscillators, etc. that are prone to 30 _ . 

misalignment, temperature caused drift, and various other blG 3 15 a ™pufied representation of a diagnostic 

modes of failure. The physical medium itself also introduces processor/mo tutor in accordance with one embodiment of 

impairments in the form of added noise and reflections. me P resent invention. 

Digital receivers can inherently correct for errors up to a FIG - 4Ade P ict s hardware for implementing a communi- 

certain threshold, but if that threshold is exceeded, commu- 3S cations unk diagnostic system in accordance with a pre- 

nication is severely impaired. ferred embodiment of the present invention. 

Digital communication techniques found their earliest FIG ' 4 . B dc P icts a s y stem architecture block diagram of a 

application in the context of point-to-point links as would be workstation. 

used by the military or commercial common carriers. In 5A * a flowchart describing the steps of analyzing 

these applications, the link may be tested by interrupting 40 a channcl response in accordance with one embodiment of 

normal service and transmitting test signals. For example, if mc P rcseQt invention. 

it is desired to test the link margin of a link, a test signal ^ IG - 5B shows a system for implicitly calculating the 

could be transmitted at progressively reduced power levels pseudo-inverse of an equalizer response in accordance with 

until the signal can no longer be accurately received. one embodiment of the present invention. 

Increasingly, digital techniques are being applied to 45 FIG - 5C * a flowchart describing steps of analyzing a 

broadcasting. One important example is broadcasting of channel response by making measurements at multiple fre- 

digital video signals to numerous subscribers over coaxial quencies in accordance with one embodiment of the present 

cable, or over the air. Commercial digital broadcast systems invention. 

require high reliability yet these digital broadcast systems 50 F* G - 5D depicts a typical structure of a channel in the 

cannot be tested as easily as the point-to-point links. Inter- context of an HFC (hybrid fiber cable) system, 

rupting normal service for transmission and measurement of FIG- 6A is a block diagram describing various modes of 

test signals is not a commercially viable option. taking advantage of error signal 217 which represents 

ingress interference in communications channel 104. 

SUMMARY OF THE INVENTION ^ FIG. 6B depicts a communications channel with an addi- 

The present invention provides method and apparatus for tional ingress signal being injected by a second digital 

non-invasive testing of digital communications test systems. transmitter. 

In accordance with a first aspect of the invention, amplitude FIGS. 7A-7B are diagrams illustrating bow system mar- 
measurements are made for multiple frequencies of a multi- gin may be tested non-invasively in accordance with one 
frequency communication system, and converted to the time go em bodiment of the present invention, 
domain. An adaptive filter then seeks to match its output to FIG. 7C is a diagram illustrating the effect of broadband 
the time domain representation, thereby characterizing the noise on a particular constellation point as received, 
channel. Impedance mismatches may be precisely located FIG. 7D is a diagram illustrating the effect of a narrow- 
using this technique. band interferer on a particular constellation point as 

In accordance with a second aspect of the invention, an 65 received, 

error signal representing a difference between a signal FIG. 7E is a diagram illustrating the effect of residual 

transmitted through the channel and a received signal is phase noise on a particular constellation point as received. 
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FIG. 7F is a diagram illustrating the effect of residual 
amplitude modulation on a particular constellation point as 
received. 

FIG. 7G is a flowchart describing steps of analyzing an 
error signal to estimate system margin and obtain other 5 
system information in accordance with one embodiment of 
the present invention. 

FIG. 7H is a diagram depicting narrowband interference 
superimposed on a digital signal. 

FIG. 71 is a diagram depicting the response of a notch 
filter for removing narrowband interference in accordance 
with one embodiment of the present invention. 

FIG. 7J is a diagram depicting the geometry of a particular 
constellation point. 15 

FIG. 7K is a diagram depicting a constellation point as 
affected by more than one type of impairment. 

FIGS. 8A-8C are diagrams illustrating how parameters 
relating to residual modulation and symbol timing variation 
may be measured in accordance with the invention. 20 

DESCRIPTION OF SPECIFIC EMBODIMENTS 
Digital Communications Background 

FIG. 1 depicts a simplified representation of a represen- 
tative digital communication system 100 suitable for testing 25 
in accordance with the present invention. Digital commu- 
nication system 100 includes a transmitter 102, a commu- 
nications channel 104, and a digital receiver 106. Transmit- 
ter 102 receives digital data to be transmitted to receiver 106 
and incorporates channel coding, modulation, 30 
up-conversion, amplification, and output filtering functions. 
The present invention is not restricted in application to any 
particular channel, code, digital modulation scheme, or 
frequency band. Channel 104 represents any physical 
medium for transferring a signal for transmitter 102 to 35 
receiver 106. For example, transmitter channel 104 may 
represent a coaxial cable delivery system and transmitter 
102 may be a cable head-end. Alternatively, channel 104 
may represent transmission between a satellite and a ground- 
based receiver. Furthermore, receiver 106 need not be 40 
uniquely coupled to channel 104. Channel 104 may repre- 
sent a broadcast medium accessible to numerous receivers. 

The present invention takes advantage of various signals 
and data generated internally to digital receiver 106. FIG. 2A 
is a block diagram depicting one implementation of a 45 
portion of digital receiver 106. After filtering and 
downconversion, a signal received from channel 104 is 
passed through a variable gain stage 202. The amplified 
output of gain stage 202 passes through an AGC control 
circuit 204 which controls the gain of gain stage 202 using 50 
a gain control signal 205. A bandpass filter 206 selects a 
particular range of frequencies so as to guard against later 
development of unwanted image signals by a tuning/ 
downconversion stage 208. Tuning/downconversion stage 
208 selects a particular signal of interest and converts this 55 
signal to baseband. Typically, tuning/downconversion stage 
208 incorporates at least one voltage controlled oscillator 
(VCO) which is at least partially controlled by a special 
control signal developed to minimize residual frequency 
modulation in the downconverted signal. For modulation 60 
schemes that involve two orthogonal components such as 
QPSK, QAM etc., tuning/downconversion stage 208 devel- 
ops two outputs corresponding to each orthogonal compo- 
nent. The present invention is not limited to any particular 
modulation scheme. 65 

Up until this stage, the signals are presumed to be analog 
and continuous. A sampling stage 210 converts the down- 
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converted signal to a discrete-time digital representation. 
The sampling rate is synchronized to symbol boundaries 
within the received data stream by a baud timing recovery 
loop 212 that generates a sampling rate control signal 213. 

An adaptive equalizer 214 automatically corrects for 
distortions in channel 104 and typically includes a digital 
FIR filter and/or IIR filter with variable tap weights. 
(Alternative receiver architectures employ equalizers with 
fixed tap weights or omit equalization entirely.) A symbol 
decision stage 216 examines the data output from adaptive 
equalizer 214 and estimates the transmitted data based on a 
maximum likelihood decision rule corresponding to the 
modulation scheme employed. A channel decoding stage 
218 removes effects of forward error correction (FEC) or 
other channel coding schemes employed. Further stages (not 
shown) may remove source coding, differential coding, 
scrambling, interleaving, or compression applied to the 
transmitted data. 

Symbol decision stage 216 also generates an error signal 
217 corresponding to the difference between an ideal modu- 
lation signal representing the output estimates and input of 
symbol decision stage 216. This error signal is an input to a 
weight updating stage 220 that adjusts weights 221 of 
adaptive equalizer 214 to minimize certain components of 
the error. A carrier recovery stage 222 also makes use of this 
error signal to remove residual frequency offset and fre- 
quency modulation effects induced in part by tuning/ 
downconversion stage 208 by controlling the internal VCO 
with a desp inning control signal 223. 

The operation of typical digital communication system 
100 will be made clearer by reference to FIG. 2B which 
depicts a so-called signal constellation for a representative 
digital modulation scheme 64 point QAM. A signal constel- 
lation is a graphical representation of the possible symbols 
for a given modulation scheme. The horizontal and vertical 
axes correspond to the orthogonal components of the modu- 
lation signal. Each possible symbol is represented by a point 
at the position of its associated (I,Q) coordinates. As shown 
in FIG. 2B, 64 point QAM is represented as an array of 64 
points. Since log2(64)-6, the choice of one particular sym- 
bol for transmission during a given symbol period commu- 
nicates 6 bits of information. 

Channel 104 introduces noise causing the constellation 
points to disperse from their ideal locations marked in FIG. 
2B. Transmitter filtering and echoing effects within channel 
104 introduce inter-symbol interference (ISI) which must be 
removed by equalization. Also, tuner/downconversion stage 
208 introduces residual frequency offset which causes the 
received signals to "spin" about the origin. 

Symbol decision stage 216 essentially receives from 
adaptive equalizer 214 a series of points in the coordinate 
space represented by the graph of FIG. 28 and assigns to 
each point the nearest constellation point as the symbol 
received. Weight updating stage 220 adjusts equalizer 214 to 
minimize the distance between the received signal and the 
constellation point transmitted. Carrier recovery stage 222 
operates to minimize the spinning of received points around 
the origin. 

FIG. 2C depicts an alternative implementation of a por- 
tion of digital receiver 106 as described in U.S. Pat. No. 
5,495,203 issued to Harp et al., EFFICIENT QAM 
EQUALIZER/DEMODULATOR WITH NON-INTEGER 
SAMPLING, assigned to the assignee of this application, 
and the contents of which are herein incorporated by refer- 
ence for all purposes. This receiver architecture will now be 
described in brief. An A-D converter 224 corresponds to 
sampling stage 210 of FIG. 2 A but does not operate at 
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baseband but rather accepts an input centered at an IF and From despinning control signal 223, equalizer weights 221, 
samples at a fixed rate determined by a reference oscillator and error signal 217, diagnostic processor/monitor 300 
226 in conjunction with a divide by 4 circuit 228. The fixed generates information about phase noise, residual frequency 
sampling rate and IF are chosen in relation to the spectrum modulation, frequency error, and frequency stability. From 
of the modulation signal so as to enable digital quadrature 5 sampling rate control signal 213 and equalizer weights 221, 
direct conversion to baseband by a digital QDC stage 230. diagnostic processor/monitor 300 generates information 
A digital rcsampler 232 synchronizes to the received data in about transmitter baud timing. Equalizer weights 221 may 
conjunction with a symbol timing phase lock loop 234 ^ e ^ to g enerate information about the channel response 
analogous to baud timing recovery loop 212. An adaptive including the time domain locations of impedance mis- 
equalizer 236 is specially adapted to operate at a sampling 10 match es- Error signal 217 generated by symbol decision 
rate that is a non-integer fraction of the symbol rate. Acarrier ^g? 1 * ?™ be ^ * 

recovery stage 238 performs the despinning function fS^S^^ " " 

directly on the digital data output from equalizer 23t iwith- FIG ' 4A d ^ faardware fof ^ lementin a communi . 

out adjustoent of oscillators as in the receiver of FIG 2A. cations liflk & sti c system 400 in accordance with a 

A symbol decision stage 240 estimates the received data is p re f err ed embodiment of the present invention. Communi- 

based on the output of earner recovery stage 238 and cations hok diagnostic system 400 includes an RF tuner/ 

generates error signals to drive earner recovery stage 238 downconverter 402, a Model 194 Wideband Signal Acqui- 

and an equalizer weight updating stage 242. sition and Playback Unit 404 available from the assignee of 

It can be seen from FIG. 2A and 2C, that digital receivers this application, and a workstation 406. Tuner/ 

generate various signals and data that correct for impair- 20 downconverter 402 receives a signal from channel 104 and 

ments in overall communication system operation. For is particularized for the band and transmission medium of 

example, equalizers 214 and 236 correct for echoing and interest Tuner/downconverter 402 selects a signal of interest 

frequency response distortions in channel 104. AGC loop and downconverts it to a predetermined IF. In the preferred 

204 inherently corrects for hum and residual amplitude embodiment, the IF is 43.75 MHz. 

modulation. The operation of baud timing recovery loop 25 Wideband Signal Acquisition and Playback Unit 404 

212, symbol timing PLL 234, and carrier recovery stages emulates the operation of the receiver components depicted 

222 and 238, inherently correct for residual frequency in FIG. 2A. The IF signal is digitized at up to 250 MHz and 

modulation, phase noise, drift, and instability, whether these captured in an internal snapshot memory. The remaining 

artifacts are introduced at the transmitter or receiver end. Of functions of the receiver of FIG. 2A are performed by 

course, if these various impairments exceed a given level, 30 workstation 406. 

the receiver will no longer be able to recover the transmitted Workstation 406 includes a processing unit 408 contain- 

data at an acceptable bit error rate (BER). ing control elements. A keyboard 410 and a mouse 412 are 

In general FIGS. 2A and 2C are merely representative of coupled to processing unit 408 and enable the system user to 

possible digital receiver architectures to which the present input commands and data to workstation 406. A display 

invention is applicable. Note that some components of 35 device 413 is provided for displaying test results. Worksta- 

digital receiver 106 are not present in every implementation. tion 406 further includes a floppy disk drive 414 adapted to 

For example, some implementations use equalizers that do receive a floppy disk 416, and an internal hard drive 418. 

not require weight updating. Other implementations do not FIG. 4B depicts a system architecture block diagram of 

employ any kind of equalizer. Many links do not employ workstation 406 of FIG. 4A. The system elements are 

channel coding to improve performance and thus do not 40 interconnected by a system bus 420 and include a display 

require a channel decoding stage. adapter 422 coupled to display device 413, an I/O controller 

Overview of Diagnostic System 424, system memory 425, a central processor (CPU) 426, a 

In accordance with one embodiment of the present speaker 428, a serial port 429, keyboard 410, a hard disk 

invention, impairments are non-invasively measured and controller 432 coupled to hard disk drive 418, and a floppy 

otherwise characterized by monitoring and analyzing inter- 45 disk controller 434 coupled to floppy disk drive 434. CPU 

mediate data and signals generated internal to a digital 426 controls operation of workstation 406. Program and data 

receiver. The interpretation and availability of these signals information are stored in system memory 424, on floppy 

will of course vary depending on the internal architecture of disk 416, and on hard disk drive 418. 

the particular digital receiver employed. Analysis software stored in system memory 424, on 

FIG. 3 is a simplified representation of a diagnostic 50 floppy disk 416, and/or on hard disk drive 418 implements 

processor/monitor 300 in accordance with one embodiment the functionality of diagnostic processor/monitor 300 as well 

of the present invention. Possible inputs to diagnostic as the functionality of digital receiver 106. Of course any 

processor/monitor 300 include gain control signal 205 from tangible medium including CD-ROM, magnetic tape, flash 

AGC control circuit 204, a despinning control signal 223 memory etc. could store this software, 

from carrier recovery stage 222, a sampling rate control 55 The workstation of FIGS. 4A-4B is merely representative 

signal 213 from baud rate recovery loop 212, the-current of hardware usable to implement the present invention. In an 

equalizer weights from adaptive equalizer 214, and the error alternative embodiment, the present invention is imple- 

signals and received symbols generated by symbol decision mented as a battery-operated portable system suitable for 

stage 216. The monitor component of diagnostic processor/ field work. 

monitor 300 receives all these signals and converts them to 60 In the embodiment shown in FIGS. 4A-4B, Wideband 

a form suitable for analysis. The inputs of FIG. 3 are Signal Acquisition and Playback Unit 404 captures sufficient 

described in reference to the receiver architecture of FIG. 2A signal information for a full analysis in less than one second, 

but similar inputs could be obtained from the receiver The analysis itself completes in 1-2 minutes. An optional 

architecture of FIG. 2C. coprocessor card installable within workstation 406 can 

From gain control signal 205, equalizer weights 221, and 65 reduce analysis time to approximately 10 seconds, 

error signal 217, diagnostic processor/monitor 300 generates In a preferred embodiment, workstation 406 derives and 

information about residual amplitude modulation and bum. displays on display device 408 combinations of the follow- 
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ing: (1) transmission (or propagation) channel response; (2) In one embodiment, the pseudoinverse is computed at 
carrier-to-noise ratio (CNR); (3) spectrum of the underlying step 502 explicitly, using Singular Value Decomposition as 
noise and any interference signals; (4) phase noise estimates; is well known to those of skill in the art. Any of a variety of 
(5) signal spectrum, signal polar plot, signal constellation well-known matrix algebraic methodologies could be used, 
and eye diagram (for analyzing intersymbol interference) 5 FIG. 5B shows a system 550 for implicitly calculating the 

and (6) other diagnostic displays. pseudoinverse at step 502. An adaptive filter 552 (distinct 

Analysis of Single Channel Response from adaptive equalizer 214) receives an ideal modulation 

In accordance with one embodiment of the present signal representing the symbol decisions generated by sym- 
invention, the response of channel 104 or a portion of bol decision stage 216. The output of adaptive filter 552 is 

channel 104 may be determined. From the channel response, 10 subtracted from the input data to adaptive equalizer 214 

the operator can determine undesirable multipath, cable, or within a summer 556. A weight updating stage 558 then 

waveguide mismatching. When channel 104 is the ambient operates to modify the tap weights of adaptive filter 552 to 

atmosphere, structures such as buildings, vehicles, aircraft, minimize the output of summer 556, preferably using the 

hills, etc., can reflect the transmitted signal such that the well-known least mean square (LMS) procedure. Since the 

transmitted signal is subject to multipath. If channel 104 is 15 filter taps have converged, i.e., adapted to a solution such 

a coaxial distribution system, reflections may be caused by that the error cannot be further minimized and the error is 

impedance mismatches resulting from loose connectors, etc. acceptably small enough, the response of adaptive filter 552 

Receiver 106 receives multiple versions of the transmitted matches the overall system response. It has been found that 

signal at slightly different times due to the differences the channel can be modeled with residual errors 35-40 dB 

between the path lengths of the reflected signals and the path 20 below the level of the channel response, 

length of the transmitted signal. Discovering such reflections Once a pseudoinverse is computed as an estimate of the 

is an important aspect of diagnosing a link. overall system response H, particular components of the 

The channel response analysis technique of the present system response may be isolated, at step 504. For many 

invention will be made clearer by a simplified mathematical communications systems, the response of the transmitter 

characterization of digital communications system 100. The 25 102, the pre-equalizer components of receiver 106, or both 

notation used herein follows a standard linear time-invariant are well known. In these cases, the known responses can be 

model. The conversion of the results to slowly time-varying deconvolved from ft, leaving an estimate of the channel 

or somewhat non-linear systems is well-known in the art. response alone. 

Referring again to FIG. 1, transmitter 102 and associated Furthermore, there may also be known portions of the 
electronics shape and transmits an input symbol stream S n . 30 response of channel 104 that may be deconvolved to focus 
Channel 104 carries the output pulses corresponding to the on a particular portion of interest. For example, it may be 
symbol stream, modifying the pulses according to the desirable to analyze only the portion of a coaxial cable 
response of the channel. Referring to FIG. 2A, the compo- system lying within a subscriber's home. Once a desired 
nents of receiver 106 to the left of adaptive equalizer 214 response is obtained, it may be shown on display device 508. 
introduce a further response component. The combined 35 In accordance with the invention, an all-pass portion of 
response of transmitter 102, channel 104, and components of the response of a channel component of interest may be 
receiver 106 to the left of adaptive equalizer 214 can be isolated at step 506 to identify the locations of reflections, 
denoted as H, while the response of adaptive equalizer 214 The time domain representation of the all-pass component 
can be denoted as W. The overall response of the system up will show the locations of individual reflections which could 
to the output of adaptive equalizer 214 can then be written 40 be caused by e.g., loose connectors. In accordance with the 
as invention, the allpass portion of the response may be 
w h s •s obtained without exact knowledge of the frequency- 
" " ^ dependent amplitude response by using homomorphic signal 
where §„+c is the output of the equalizer structure and is a processing or parametric signal processing, 
noisy estimate of the original symbol stream §„+€. 45 Homomorphic signal processing is a method by which a 
When communications system 100 is working effectively, signal is decomposed into an all-pass component and a 
S„+e is a very good estimate of S m i.e., the received symbols minimum-phase component. This decomposition can be 
are very close to the same symbols transmitted. In this case, performed using techniques well known to those skilled in 
W operating on H as a matrix operator is very close to the me art » including methods described in [Oppenheim89] and 
identity, I, a matrix operator being defined as an operator 50 [IEEE79], the contents of which are herein expressly incor- 
which takes one or more vectors, matrices, or tensors as porated by reference for all purposes. In the case where the 
arguments and returns a vector, matrix, or tensor. Thus, the input signal consists of a signal corrupted by micro- 
following relationships hold: reflections, the minimum-phase portion of the signal will 

contain a representation of the spectrum of the direct arrival 

S """^» ( 2 ) 55 of the desired signal, while the all-pass portion will contain 

w H s„-s„±5W h*l (3) a representation of the impulse response of the charmers 

micro-reflections. 

It then follows that the pseudoinverse of W, W#, is very fa cases where the ideal (uncorrupted) channel response is 

close to H, the overall system response. One can estimate H known, or can be reasonably estimated, parametric signal 

as: 60 processing techniques may be employed to better estimate 

H„w#. ^ the corrupted channel response. In using parametric 

approaches, a model for the corrupted channel response is 

FIG. 5A is a flowchart describing the steps of analyzing hypothesized. This model includes several unknown 

a channel response in accordance with one embodiment of parameters, such as time-of-arrival and amplitude of micro- 

the present invention. At step 502, the pseudoinverse of the 65 reflections. The values for each of these parameters may 

equalizer response is determined using the weights gener- then be determined in any number of methods well known 

ated by weight updating stage 220. to those skilled in the art, including Maximum Likelihood, 
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Maximum Entropy, Minimum Variance, etc. Details of these available. This defines the response of the equalizer with a 

methods, as known to those of skill in the art, may be found resolution of T/6 (<0.17 T); three times better than that 

throughout the literature, including [VanTrees68], obtainable with 172 equalizers and their inherent 0.5 T 

[Haykin83], [Haykin85], [Haykin91], and [Marple87], all resolution. 

the contents of these references being herein expressly 5 Analysis of Multi-Frequency Multi-Channel Response 

incorporated by reference for all purposes. Even fi°« resolution in identifying reflection sources is 

Thus, the locations of individual reflections may be obtainable by making measurements at multiple frequencies 

obtained without resort to invasive probing of the system across me bandwidth of the entire cable system as opposed 

while it is operating. This is an important advantage over 10 monitorin g a signal. The cable system bandwidth is 

prior art time domain reflectometry techniques which 10 200 MHz tc \ 1 F ° r a 5( ^. *JHz s y stem > 

require that normal operation be interrupted. reflection sources may be localized to within 0.7 m. 

THe resolution of the above^iescribed system in locating h mQ \ 5C 15 a fl ° Wch ^ desCTlbm S ° f ™£yf* a 
« • t« •« j u *u «? *• * i channel response by making measurements at multiple fre- 

reflection sources is limited by the effective temporal reso- quencies inVccordance with one embodiment of the present 
ktion of the equalizer For a simple finite impulse response Ai step 560 , the preferred embodiment measures 

(FIR) structure, the effective temporal resolution is directly 15 envelope amplitude for each signal carried by the cable 
related to the tap spacing. If standard processing approaches system. Typically, a measurement is made every 6 MHz by 
are utilized, a typical T- or T/2-spaced equalizer is limited to tuning receiver 106 to each channel. The amplitude mea- 
a resolution of T or T/2, respectively, where T is equal to the surement must be made at a stage of receiver 106 prior to the 
symbol period. For typical digital cable television (CATV) circuitry shown in FIGS. 2Aand 2C, e.g., at the input to gain 
applications, a 5 MBaud, T/2-spaced equalizer would have 20 stage 202. The result of step 560 is a list of frequencies and 
a temporal resolution of 100 ns. This translates to a resolu- measured amplitudes. 

tion of about 24 meters (m) in typical coaxial cable. At step 562, the list of frequencies and measured ampli- 

In accordance with one embodiment of the invention, by tudes is converted into the time domain by applying an FFT 
employing a non-integer, fractionally-spaced equalizer, the (Fast Fourier Transform). This results in a time domain 
temporal resolution can be increased by interpreting the 25 representation of the signal measured over an aggregation of 
equalizer as a polyphase structure. In such systems, the the frequencies. At step 564, an impulse is repeatedly fed 
sample rate is approximately equal to (2*N*K/L) times the into an adaptive digital filter. The adaptive filter operates to 
symbol rate where 2*L>K>L, and K, L, and N are integers. modify its transfer function to force its output to equal the 
In the case of a T/1.2-spaced equalizer, as is described in the time domain representation developed in step 562. The 
Harp et al. patent, the resulting temporal resolution is T/6. 30 adaptation method is preferably the well-known gradient 
For a 5 MBaud system, this would provide a 33 ns descent algorithm. In one embodiment, the adaptive digital 
resolution, better than 8 m in the same coaxial cable where filter is an all-pole filter, and adaptation includes positioning 
only 24 m was the resolution available using a T/2-spaced only the poles of the transfer function. In another 
equalizer. Other spacings can provide even better resolution embodiment, the adaptive digital filter includes zeroes as 
than that of a T/1.2 equalizer, though usually at the expense 35 well which are positioned during adaptation, 
of greater computational complexity. The adaptation process forces the response of the adaptive 

With a T/1.2-spaced equalizer, one can visualize the digital filter to approximate the response of channel 104 as 
structure as a linear, time-varying system comprised of five measured across the entire system bandwidth. This is an 
distinct equalizer phases. Each of these phases is used in turn example of parametric signal processing. At step 566, the 
as each new sample is received. Without loss of generality, 40 adaptive digital filter response is converted to the time 
the first of these phases (Phase 0) is assigned a relative delay domain. If the filter response includes both poles and zeroes, 
of 0 and is responsible for providing an output correspond- only the "poles" component of the response need be con- 
ing to some symbol sent at t«0. The next phase (Phase 1) of verted to the time domain for this purpose. The resulting 
the equalizer is then responsible for predicting the symbol time domain representation of the poles or "resonances" will 
sent at t=T, but using samples collected with a relative delay 45 show pulses corresponding to the locations of impedance 
of 5T/6. Thus the taps of Phase 1 are effectively T/6 mismatches, as caused by, e.g., loose or corroded connec- 
(T-5T/6) advanced relative to those in Phase 0. Similarly, tors. 

the remaining phases are advanced by T/6 from the previous At step 568, if the filter response includes zeroes as well, 
phase. This relationship may be summarized as in Table 1: the "zeroes" component of the response is isolated and 

50 displayed in the frequency domain. This shows the effects of 
TABLE 1 shaping filters in transmitter 102, fitters used at intermediate 

points in channel 104 and any other filter frequency response 
distortions along the way. 
Localization of Channel Impairments 
55 Generally, accurate identification of the location of an 
impedance mismatch requires isolation of a portion of the 
response of channel 104. FIG. 5D depicts a typical structure 
of channel 104 in the context of an HFC (hybrid fiber cable) 
system. A fiber-optic transmitter 570 transmits optical 
60 energy to an electro-optical node 572 Electro-optical node 
572 transmits signals through a distribution system 576 to a 
tap box 578. At tap box 578, an electrical signal is tapped off 
As shown in this table, the relative delay decreases by T/6 for transmission through drop wiring 579 to premises 580 of 
from phase to phase, and increases by 5T/6 from tap to tap. a particular subscriber. Premises 580 include televisions 582 
As can also be seen from this table, by properly ordering the 65 and a set-top unit 584 implementing receiver 106. 
N taps of the entire equalizer structure, all relative delays In FIG. 5D, nodes A-E represent points where the signal 

from -4T/6 up through (N-5) T/6 in increments of T/6 are may be tapped off and monitored to identify the channel 
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response up to that point, by any of the procedures discussed or equalizer 236) and is generally considered to be the power 

in reference to FIGS. 5A-5C. To isolate a particular segment of the signal divided by the sum of the power of the in-phase 

of channel 104, responses are obtained at a node following error plus the power of the quadrature error. Alternatively, 

the segment and a node preceding the segment. To obtain the the cluster variance may be measured based on a modulated 

response of the optical system represented by fiber-optic 5 si S nal carrying a reencoded version of the data output from 

transmitter 570 and electro-optical node 572, the response channel decoding stage 218. The cluster variance represents 

obtained at node A is deconvolved from the response an estimate of the signal-to-distortion power ratio. Depend- 

obtained at node B. To obtain the response of distribution 5« on the Portion type and severity, the bit error rate 

system 576, the response obtained at node B is deconvolved £* ER ) can be calculated as a function of the signal-to- 

from the response obtained at node C To obtain the response 10 distomon power raUo or chister vanance. 

r . „ • C7Q *ho «,™„c-o ^M.-no^ ot n ^/n .v Precise calculation of the signal margin (the difference 

of drop wiring 579 the response obtained at node C , s between ^ si (S ^ } ^ a Tthe required BER 

deconvolved from the response obtained at node D. To ys ^ for a given impairment) in accordance 

obtain the response of the wmng inside subscriber premises, ^ the m mvenlion ^ disced be low. However, the 

the response obtained at node D is deconvolved from the prcscnt invenlion ^ pmvidcs for a simplified calculation 

response obtained at node E. is Q f BER and signal margin directly from the cluster variance 

Analysis of Error Signal based on the recognition that in most cases distortion can be 

In accordance with the invention, diagnostic processor/ considered to have an approximately Gaussian distribution 

monitor 300 receives error signal 217 corresponding to the and be independent. The Gaussian assumption holds gener- 

difference between an ideal modulation signal representing ally in digital cable systems since the equalizers and other 

the output estimates and the input of symbol decision stage 20 filtering elements operate as "summers" and in accordance 

216. For systems that employ channel coding, an alternative with the well-known Central Limit Theorem cause arbitrary 

way to generate the error signal is to reencode the decoded probability distributions to become Gaussian. The indepen- 

output of channel decoding stage 218 and form an ideal dence assumption also holds for digital cable systems 

modulation signal representing this corrected encoded data. because they generally incorporate large interleavers or error 

The error signal is then the difference between this encoded 25 randomizers that break up burst errors (or distortions that 

data and the input of symbol decision stage 216. might cause burst or cause dependence of the errors) and 

FIG. 6A is a block diagram describing various modes of randomize them resulting in independence of error events, 
taking advantage of error signal 217 which represents Bit Error Rate Based on Error Signal Analysis 
ingress interference in communications channel 104. One One embodiment of the present invention makes use of 
possibility is direct display of the error signal in the time 30 the above two assumptions and approximates the signal-to 
domain. An operator able to view this display may be able noise ratio (SNR) to be the cluster variance. Each modula- 
te recognize the source of ingress error as for example, 60 rion scheme has an associated BER vs SNR graph which is 
Hz line noise. Burst errors or other time varying phenomena used to estimate the BER from the measured cluster van- 
may be identified. Alternatively, it may be useful to convert ance. The signal margin may be estimated similarly. For 
the signal into the frequency domain (block 602) using 35 systems that make use of error correction, the post-decoded 
well-known techniques prior to display. This provides the bit error rate (DBER) is estimated by first computing the 
equivalent of a spectrum analyzer on an unoccupied channel. pre -decoded error rate from the cluster variance and apply- 
A frequency domain display may show the presence of ing the input vs output error rate curves for a given decoder 
ingress interference from an AM or FM broadcast signal to compute the output error rate of the decoder. Curves for 
leaking into the cable system, for example. 40 a variety of error correcting codes are efficiently stored in a 

The ingress display can also be normalized to the signal look-up table within block 604. There is no need to simulate 

power to display the power spectral density of the ingress the operation of the decoder. 

signal relative to the signal power. The signal power may be Coarse System Margin Approximation Based on Error Sig- 

obtained at the output of adaptive equalizer 214 or 236. This nal Analysis 

signal power then is used to generate a reference power level 45 Avoiding simulation of decoder operation produces great 

on the frequency display (e.g., 0 dBc or decibels relative to benefits in determining margin under low error rate condi- 

the signal or carrier power). The ingress power is then tions. This is because measuring low error rates can be very 

plotted (block 604) as a spectral density relative to the signal long and tedious. With powerful error-correction coding, 

power (e.g., -40 dBc/Hz) or as a power in the resolution reasonably short measurements will produce either zero 

bandwidth of the display. For example, if the resolution 50 errors above the SNR threshold or a 50% bit error rate below 

bandwidth is 1 kHz, the -40 dBc/Hz above is multiplied by the SNR threshold, due to the breakdown of the error- 

1000 Hz (30 dB) to obtain the ingress power of -10 dBc. correction coding scheme. This gives very little indication of 

The ingress display can be further corrected for adjust- available margin. However, subtracting a required system 

ments made in the equalization by adjusting in the frequency SNR from the SNR estimate derived from the cluster 

domain for the influence of the equalizer spectral response 55 variance yields a useful estimate of the available margin. A 

on the ingress signals to determine the original ingress technique for obtaining a more accurate estimate of system 

power prior to equalization processing. margin is discussed below. 

Alternatively, statistics of error signal 217 may be com- Precise Evaluation of System Margin by Analysis of Error 

puted (block 604) in relation to the symbol stream output Signal 

from symbol decision stage 216 to develop indicia of link 60 A more refined estimate of system margin may be 

quality. One extremely useful measure of quality is the obtained in accordance with one embodiment of the present 

"cluster variance" (CV). The cluster variance is the ratio of invention. This is done by separating the component of the 

the variance of the power in the error signal to the power in error signal which is caused by broadband noise from the 

the symbol stream. components caused by narrow-band interferers, residual 

Cluster Variance Based on Error Signal Analysis 65 amplitude modulation, and residual phase modulation. 

The cluster variance is measured around the demodulated FIG. 7C depicts a constellation point 756 as affected by 

signal constellation point (e.g., at the output of equalizer 214 broadband noise only. Ideally, constellation point 756 would 



07/08/2004, EAST Version: 1.4.1 



US 6,385,237 Bl 



13 



14 



be a dimensionless point but the effect of broadband noise is 
to diffuse the received constellation point into a disk. FIG. 
7D depicts constellation point 756 as affected by a narrow- 
band interferer. The effect of the narrowband interferer is to 
create a ring 758 having the ideal location of constellation 5 
point 756 as a center. 

FIG. 7E depicts constellation point 756 as affected by 
residual phase modulation. FIG. 7E assumes that constella- 
tion point 758 is in either the upper right or lower left 
quadrant of a constellation such as the constellation in FIG. 10 
2B and located on a line through the origin at roughly a 45 
degree angle to the axes. The effect of residual phase 
modulation is to disperse the received constellation point in 
a direction tangent to a circle including constellation point 
756 and having the origin of FIG. 2B as a center. It should 15 
be noted that the residual phase modulation of FIG. 7E will 
not include errors corrected for by carrier recovery stages 
222 or 238. 

FIG. 7F depicts constellation point 756 as affected by 
residual amplitude modulation. The same constellation 20 
geometry assumptions made in reference to FIG. 7E also 
apply to FIG. 7F. The effect of residual amplitude modula- 
tion is to disperse the received constellation point in a 
direction perpendicular to a circle including constellation 
point 756 and having the origin of FIG. 2B as a center. It 25 
should be noted that the residual amplitude modulation of 
FIG. 7F will not include the errors corrected for by AGC 
stage 204. 

For an actual received signal, the errors of FIGS. 7C-7F 
will be combined. The present invention provides a tech- 30 
nique for isolating these errors. 

FIG. 7G is a flowchart describing steps of analyzing an 
error signal to estimate system margin and obtain other 
system information in accordance with one embodiment of 
the present invention. At step 760, the preferred embodiment 35 
analyzes the spectrum of the received signal as would be 
present, e.g., on the output of carrier recovery stage 238 or 
adaptive equalizer 214 to identify interferers. FIG. 7H 
depicts such a spectrum 762 showing two narrow band 
interferers 764 and 766. Any narrowband signal having 40 
greater than three times the median magnitude of spectrum 
762 is considered to be an interferer. At step 768, a notch 
filter as shown in FIG. 71 is applied to remove the narrow 
band interference. Referring to FIG. 71, a notch filter 
response 770 is constructed so as to remove the interferers 45 
identified in step 760. Notch filter response 770 includes 
nulls 772 and 774 corresponding to the frequency and 
bandwidth of narrowband interferers 764 and 766 respec- 
tively. FIG. 71 of course shows ideal filter performance. 

At step 772, the preferred embodiment performs statisti- 50 
cal analysis on samples from the two-dimensional error 
signal 217 corresponding to a selected constellation point, 
an in-phase error correlation signal E {c/€ 7 } corresponding 
to variance along the in-phase (horizontal) axis, a quadrature 
error correlation signal E {e^Co} corresponding to variance 55 
along the quadrature (vertical) axis, and an error cross- 
correlation signal E {e,e G } corresponding to the cross- 
correlation between error samples along the in-phase and 
quadrature axes. 

At step 774, the preferred embodiment determines a 2 AAf , 60 
corresponding to the amount of residual AM, corre- 
sponding to the amount of residual PM, and o 2 ^ correspond- 
ing to the amount of wideband noise based on the correlation 
signals computed in step 772. Step 774 takes advantage of 
a relationship among E^,}, E^e^}, E{eje Q }, o 2 ^, 65 
cPpM, and <r N determined by the position of the selected 
constellation point in relationship to an origin 776. Referring 



now to FIG. 7J, constellation point 756 lies at a distance r 
from the origin and a vector 778 from origin 776 to 
constellation point 756 is at an angle 6 to a horizontal axis 
780. Accordingly, the following relationships will hold: 



E{£^ Q }=r 2 sinecos9(o 2 ^ A/ -o 2 p, M )where 



(1) 
(2) 
(3) 



o 2 , and o 2 ^ are the in-phase and quadrature components of 
the broadband noise respectively, and are each equal to o 2 ^ 
The preferred embodiment solves for o 2 ^, 0%^, and 
using the above equations. 

At step 782, the preferred embodiment uses its knowledge 
°f ^am> ^pm* aQ d cr 2 ^ to compute the bit error rate. 
Referring now to FIG. 7K, a decision region 784 surround- 
ing ideal constellation point 756 is shown. Constellation 
point 756 is shown with a combination of mostly residual 
phase modulation with some broadband noise added. As 
long as the constellation point representing the received 
noisy signal falls within the boundaries of decision region 
784, the originally transmitted symbol will be correctly 
identified by the receiver. To evaluate the probability of an 
error occurring, the preferred embodiment runs a Monte 
Carlo simulation where the ideal transmitted constellation 
point is summed with pseudo-randomly chosen values cor- 
responding to simulated broadband noise, residual ampli- 
tude modulation, and residual frequency modulation with 
distributions defined in part by the o 2 ^^ and o 2 ^ 

values identified in step 780. The broadband noise and 
residual phase modulations are assumed to have Gaussian 
distributions variances cr 2 ^, and o 2 ^ respectively. The 
residual amplitude modulation is modeled as an (e.g., 120 
Hz) periodic random process with an amplitude defined by 
^am- A simulated received point falling outside decision 
region 784 is considered an error. The bit error rate is 
determined in this way and displayed. The bit error rate 
simulation may be refined by taking into account the effects 
of channel coding by applying simulated channel coding and 
decoding. 

At step 786, the preferred embodiment evaluates the 
system margin, determining the amount of extra broadband 
noise, residual amplitude modulation, and residual phase 
modulation that can be tolerated before the bit error rate rises 
above an unacceptable threshold. This is done by repeating 
the Monte Carlo simulation described in reference to step 
782 for combinations of elevated values of o 2 ^, a 2 ^, and 
o 2 ^. The combinations for which the bit error rate is 
equivalent to the threshold define a three-dimensional sur- 
face which may also be displayed to graphically illustrate to 
the operator bow much further system degradation is toler- 
able. 

In summary, the system margin is calculated as the 
amount of noise power that can be added while maintaining 
BER at its threshold value while considering the cumulative 
effect of other impairments. The present invention, however, 
also provides estimates of the degree to which non-noise 
impairments may be increased before bringing BER above 
threshold. 

Demodulation of the error signal (block 606) may also be 
useful. Converting the demodulated signal to audio may 
allow ready identification of broadcast interference. 
Alternatively, the demodulated error signal may be used to 
recover one or more deliberately injected in-band ingress 
carriers. If the level of these carriers is sufficiently low (e.g., 
-27 dBc for 64-QAM), the digital carrier is not appreciably 
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disturbed. FIG. 6B depicts communications channel 104 
with an additional ingress signal being injected by a second 
digital transmitter 620. 

Evaluation of System Margin by Applying a Synthetic 
Impairment 5 

FIGS. 7A-7B are diagrams illustrating how system mar- 
gin may be tested non-invasively in accordance with one 
embodiment of the present invention. FIG. 7 A depicts one 
approach to testing system margin. A summer 702 is inserted 
prior to adaptive equalizer 214 to add in an impairment to 
signal. The output of symbol decision stage 216 can then be 
analyzed to measure the deterioration in performance that 
results from the impairment signal in the manner described 
in connection to FIG. 6A. 

An impairment generator 704 generates the impairment 15 
from error signal 217 generated by symbol decision stage 
216. One approach is simply to amplify error signal 217 
which represents the existing interferers within impairment 
generator 704 until the bit error rate becomes unsatisfactory. 
The degree of amplification necessary to produce this level 20 
of signal quality degradation then represents how much 
more of the current interference can be tolerated. 
Alternatively, filtering could be applied within impairment 
generator 704 to isolate a particular type of interference to 
be re- applied to the signal and determine how much more of 25 
this interference could be tolerated. Another alternative is to 
apply additive noise through summer 702 and determine 
how much noise is necessary to disrupt communication. One 
or more reflections may be simulated by using impairment 
generator 704 to delay a reconstruction of the transmitted 30 
signal based on the symbol estimates. Also, if desired, only 
a particular spectral portion of the reconstruction may be 
delayed in this way. 

FIG. 7B depicts an alternative approach to estimating 
system margin. An impairment processor 750 may apply 35 
various form of signal impairment to a reconstruction of the 
transmitted signal based on the output of symbol decision 
stage 216. If forward error correction is being used, impair- 
ment processor 750 may instead use the corrected data as 
input and re-apply the relevant code prior to further pro- 40 
cessing. Impairment processor 750 may apply additive 
impairments as in FIG. 7A or may apply a simulated channel 
response. The output of impairment processor 750 is applied 
to an equalizer 752 similar to equalizer 214. (The updating 
of the weights of equalizer 752 is not shown.) Equalizer 752 45 
thus seeks to correct for the impairments introduced by 
impairment processor 750. A symbol decision stage 754 then 
makes symbol decisions based on the output of equalizer 
752. Together equalizer 752 and symbol decision stage 754 
constitute an auxiliary demodulator. The effects of the 50 
impairments introduced by impairment processor 750 are 
analyzed based on error and symbol estimate outputs of 
symbol decision stage 754 as described in reference to FIG. 
6A. 

Impairment processor 750 may apply additive interfer- 55 
ence based on error signal 217 in a manner similar to that 
discussed in reference to FIG. 7A. As the interference level 
is varied, the error and symbol estimate outputs of symbol 
decision stage 754 may be monitored for unacceptable 
performance deterioration. 60 

Impairment processor 750 may also apply a digital filter 
that simulates an intensification of an existing distortion. 
Impairment processor 750 may determine the existing dis- 
tortion by taking the pseudoinverse of the weights 221 
applied by weight updating stage 220. If a particular 65 
distortion, e.g, a reflection is to be applied, it may be isolated 
by use of simple deconvolution, homomorpbic signal 
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processing, or parametric signal processing as was discussed 
with reference to FIG. 5A. Varying the intensity level of the 
simulated reflection or reflections while monitoring the 
outputs of symbol decision stage 754 helps determine how 
much more impairment can be tolerated without disrupting 
communications. 

Analysis of Receiver Loop Parameters 

FIGS. 8A-8C are diagrams illustrating how parameters 
relating to residual modulation and symbol timing variation 
may be measured in accordance with the invention. FIG. 8 A 
depicts how parameters relating to residual amplitude modu- 
lation may be measured and displayed in accordance with 
one embodiment of the present invention. Residual ampli- 
tude modulation may be caused by either transmitter or 
receiver components. Referring to receiver 106, residual 
amplitude modulation will make itself apparent in several 
ways. One way is variation in gain control signal 205 which 
will adjust so as to remove residual amplitude modulation 
within its loop bandwidth. The variation in the overall gain 
of equalizer 214 as controlled by weight update stage 220 
will operate to remove amplitude variations seen within 
adaptive equalizer 214. Amplitude modulation that is not 
removed by the operation of AGC stage 204 and adaptive 
equalizer 214 will appear as a radial component in error 
signal 217 generated by symbol decision stage 216. 

In order to generate information about amplitude 
modulation, an amplitude modulation processing stage 802 
receives gain control signal 205, the weights generated by 
weight updated stage 220 and error signal 217 generated by 
symbol decision stage 216. Prior to further processing, gain 
control signal 205 is filtered to remove the known dynamics 
of the gain control loop. The received weights are processed 
to isolate the time varying portion of the overall gain of 
equalizer 214. The cutoff frequency of this isolated gain 
response is determined responsive to the loop bandwidth of 
the AGC loop formed by gain stage 202 and automatic gain 
control stage 204. Two-dimensional error signal 217 is 
processed to isolate a radial component. Amplitude modu- 
lation stage 802 then determines, in the time domain, a 
weighted sum of gain control signal 205, the time varying 
portion of the gain of equalizer 214, and radial component 
of error signal 217 to obtain a time domain representation of 
the residual amplitude modulation of the received signal 
input to gain stage 202. One of skill in the art will understand 
how to adjust the weights appropriate to the parameters of 
adaptive equalizer 214 and symbol decision stage 216. 

Once a time domain representation of the residual ampli- 
tude modulation is generated, it may be displayed directly. 
Alternatively, it may be converted to the frequency domain 
using well known techniques prior to display as an ampli- 
tude modulation spectrum. This type of display is highly 
useful in isolating, e.g., hum caused by poor isolation 
between power line noise and transmitter or receiver com- 
ponents. Another useful parameter that may be displayed 
numerically is the percentage of amplitude modulation 
which is defined as the amplitude of the time domain 
amplitude modulation signal divided by the overall received 
signal amplitude, which is a fixed parameter of the receiver. 

FIG. 8B depicts how parameters relating to frequency 
and/or phase modulation may be determined in accordance 
with one embodiment of the present invention. Referring 
again to receiver 106, residual frequency and/or phase 
modulation will show up in despinning control signal 223 as 
it seeks to correct for the angular component of error signal 
217. Equalizer 214 will also seek to correct for residual 
frequency and/or phase modulation and variation through 
variation in the phase component of the weights generated 
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by equalizer weight update stage 220. However, the remain- demodulator, system parameters, modulation types, displays 

ing uncompensated angular component of error signal 217 and diagnostics to meet specific system requirements, 

will also represent a component of frequency and/or phase Id another preferred embodiment, the invention is pow- 

modulation. ere d by a battery, increasing the portability of the system. 

In accordance with one embodiment of the present 5 Variations on the methods and apparatus disclosed above 
invention, a frequency modulation processing stage 804 ^ be readil y apparent to those of ordinary skill in the art. 
receives despinning control signal 223, the weights gener- Accordingly, particulars disclosed above are intended not as 
ated by weight update stage 220, and error signal 217. Prior Rations of the invention, 
to further processing, despinning control signal is filtered to ? T ,^ PJ^ _ _ „ „ „. w , , 
remove the known characteristics of the Carrier recovery 10 JJJg^l ^y^Non-L^ar Methods of Spec- 
loop Frequency modulation processing stage 804 isolates "r^kmSS] S^TS. Jk£ Array Signal Processing, 
the time varying phase component of the weights and the (Prentice-Hall 1985) 

angular component of error signal 217 and generates a [ H aykin91] S. S. Haykin, Advances in Spectrum Analysis 

weighted sum of these components with despinning control and Anay Processiri ^ prentice-Hall 1991) 

signal 223. The weighted sum serves as a time domain is [IEEE79] Programs for Digital Signal Processing, Edited 

representation of instantaneous phase error. by the Digital Signal Processing Committee, IEEE 

This phase error signal can be displayed directly in the Acoustics, Speech, and Signal Processing Society, (IEEE 

time domain, enabling the operator to recognize periodicity Press 1979). 

or other patterns. The preferred embodiment also has the [Marple87] S. L. Marple, Digital Spectral Analysis, 

capability of converting the phase error signal to the fre- 20 (Prentice-Hall 1987). 

quency domain. The carrier signal, whose amplitude and [Oppenheim89] A- V. Oppenheim & R. W. Schafer, 

frequency are fixed receiver parameters, can be superim- Discrete-Time Signal Processing, (Prentice-Hall 1989). 

posed over a frequency domain display of the phase error to [VanTrees68] E. L. Van Trees, Detection, Estimation, and 

provide convenient read out of phase noise level in dBc/Hz Linear Modulation Theory, (John Wiley & Sons 1978). 

(decibels below carrier per Hertz). The phase error signal 25 ^ e c ^ a i m: 

can also be differentiated to provide a time domain display 1 In a di S ital communication system wherein a modu- 

of frequency error relative to expected carrier frequency. late f J«°f l 3 s transmitted over a channel to a receiver, the 

Low frequency variations in frequency error represent drift mooted sipal comprising a data signal modulated onto a 

ofcarrierfrequency.Astabilitypercentageoveragiventime c ^ r a method of diagnosing said channel com- 

may also be computed and displayed. An average frequency 30 P 115 " 1 ^ f s e P s ° ■ 

error can also be displayed numerically as a fixed offset a > £^! g 531(1 moduJated S1 8 nal 35 corru P ted b y 531(1 

indicating, e.g., transmitter component variation. Note that . x ' ... . 

some receiver architectures do not downconvert received b > g^^S ™ ^ **n»te representative of an 

• « i, . w _ a~ error between said modulated signal as transmitted and 

signals all the way to DC. in these architectures, down con- . , , « . , . , 

» j n. rr . . , , , ■, . V , . , said received modulated signal: 

verted frequency onset should be removed poor to display 35 . A . . . 

rr _ J A r , _ - r*u c ) generating a power estimate representative of a power 

or frequency error. A frequency domain representation of the * -j • j j. i . j ■ i 

j.-c ^ . , ^ . . .„ . , . r of said received modulated signal; 

differentiated phase error signal will help show sources of ,v . . , . . °^ . - 

residual FM d) converUng said error signal estimate to a frequency 

FIG. 8C depicts how baud rate and baud rate variation x domain representation, and 

may be analyzed and displayed in accordance with one 40 e > Asp^g frequency domam representation rela- 

embodiment of the present invention. Variations in baud Uve t0 ^ P° wer 

time will appear within receiver 106 via sampling rate 2. The method of claim lwherem said receiver comprises 

control signal 213 which operates to synchronize sampling an *J ualizer and forther comprising the step of: 

and via the variation in delay through equalizer 214. A baud adjusting said frequency domain representation for a 

timing processing stage 806 receives sampling rate control 45 response of said equalizer. 

signal 213 and the weights generated by equalizer weight }• ^ method of claim * wherein step b) further corn- 
update stage 220. Prior to further processing, baud timing pnses 

processing stage 806 filters sampling rate control signal 213 employing symbol decision to identify an estimate of the 

to remove the known dynamics of the baud timing recovery transmitted data signal from a demodulated version of 

loop. Baud timing processing stage 806 isolates the time so me received modulated signal; 

varying delay component of the response of equalizer 214 as generating the error signal estimate using the estimate of 

represented by the weights and obtains a weighted sum of the transmitted data signal and the demodulated version 

sampling rate control signal 213 and the time domain of the received modulated signal. 

representation of the isolated delay component to form a 4. The method of claim 3 wherein step b) further 

timing signal. The so-obtained timing display may be dis- 55 comprises, prior to performing symbol decision, performing 

played numerically as the current baud rate or graphically to adaptive equalization on the demodulated version of the 

show variation in baud rate over time. Baud rate variations received modulated signal. 

would be caused by errors or problems in the transmitter. A 5. An apparatus for diagnosing a channel in a digital 
baud rate stability percentage over a given time may also be communication system, wherein a modulated signal is trans- 
displayed. 60 mitted over the channel to a receiver, the modulated signal 
Of course, in receiver architectures that do not include an comprising a data signal modulated onto a carrier signal, the 
adaptive equalizer, estimates of the above -described modu- apparatus comprising: 

lation and timing parameters could be obtained from the a) a receiver operable to receive said modulated signal as 

various control loop parameters alone. Equalizer weights corrupted by said channel, 

and error signals would not be considered. 65 generate an error signal estimate representative of an 

In a preferred embodiment, the invention is sufficiently error between said modulated signal as transmitted 

modular and programmable to allow customization of the and said received modulated signal, and 
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generate a power estimate representative of a power of 
said received modulated signal; 
b) a diagnostic processor operably connected to receive 
said error signal estimate from the receiver, the diag- 
nostic processor operable to convert said error signal 
estimate to a frequency domain representation and 
determine said frequency domain estimate relative to 
said power estimate. 

6. The apparatus of claim 5 further comprising a display 
operably connected to the diagnostic processor to receive the 
frequency domain representation relative to said power 
estimate and generate a visible display thereof. 

7. The apparatus of claim 5 wherein a processor performs 
at least some of the operations of the receiver and at least 
some of the operations of the diagnostic processor. 

8. The apparatus of claim 5 wherein the receiver includes 
a processor, said processor operable to: 

perform symbol decision on a demodulated version of the 

received modulated signal to generate an estimate of 

the data signal as transmitted; and 
estimate the error signal using the estimate of the data 

signal as transmitted and the demodulated version of 

the received modulated signal. 

9. The apparatus of claim 8 wherein the processor is 
further operable to perform adaptive equalization on the 
demodulated version of the received modulated signal prior 
to performing symbol decision on the received modulated 
signal to generate the estimate of the data signal as trans- 
mitted. 

10. The apparatus of claim 9 wherein the diagnostic 
processor is further operable to adjust said frequency 
domain representation for a response of said adaptive equal- 
ization. 

11. The apparatus of claim 5 wherein the receiver further 
includes an analog tuning/down conversion circuit and a 
sampling circuit. 

12. The apparatus of claim 5 wherein the receiver further 
includes a digital down conversion circuit. 

13. A computer program product that controls a processor 
of a receiver, the processor operable to diagnose a channel 
in a digital communication system wherein a modulated 
signal is transmitted over the channel to the receiver, the 
modulated signal comprising a data signal modulated onto a 
carrier signal, the receiver operable to receive the modulated 
signal as corrupted by said channel, the computer program 
product comprising: 

a) code to generate an error signal estimate representative 
of an error between said modulated signal as transmit- 
ted and said received modulated signal; 

b) code to generate a power estimate representative of a 
power of said received modulated signal; 

c) code to convert said error signal estimate to a frequency 
domain representation; and 
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d) code to determine said frequency domain representa- 
tion relative to said power estimate. 

14. The computer program product of claim 13 further 
comprising code to cause a display to display said frequency 

5 domain representation relative to said power estimate. 

15. The computer program product of claim 13 wherein 
said receiver includes an equalizer and the computer pro- 
gram product further comprises code to adjust said fre- 
quency domain representation for a response of said equal- 

10 izer. 

16. The computer program product of claim 13 compris- 
ing: 

code to perform symbol decision to identify an estimate of 
l5 the transmitted data signal from a demodulated version 
of the received modulated signal; and 
code to generate the error signal estimate using the 
estimate of the transmitted data signal and the demodu- 
lated version of the received modulated signal. 
20 17. The computer program product of claim 16 further 
comprising code to, prior to performing symbol decision, 
perform adaptive equalization on the demodulated version 
of the received modulated signal. 

18. An apparatus for monitoring ingress on a channel in 
25 a digital cable communication system, wherein a modulated 
signal is transmitted over the channel to a receiver, the 
modulated signal comprising a data signal modulated onto a 
carrier signal, the apparatus comprising: 

a) the receiver operable to receive said modulated signal 
30 as corrupted by ingress on said channel, 

generate an error signal estimate representative of an 
error between said modulated signal as transmitted 
and said received modulated signal; 

b) a diagnostic processor operably connected to receive 
35 said error signal estimate from the receiver, the diag- 
nostic processor operable to convert said error signal 
estimate to a frequency domain representation, the 
frequency domain representation representative of 
ingress on said channel. 

40 19. The apparatus of claim 18 further comprising a 
display operably connected to the diagnostic processor to 
receive the frequency domain representation from the diag- 
nostic processor and generate a visible display thereof. 

20. The apparatus of claim 18 wherein the receiver is 
45 further operable to generate a power estimate representative 

of a power of said received modulated signal, and wherein 
the diagnostic processor is further operable to determine said 
frequency domain estimate relative to said power estimate. 

21. The apparatus of claim 20 further comprising a 
50 display operably connected to the diagnostic processor to 

receive the frequency domain representation relative to said 
power estimate and generate a visible display thereof. 

***** 
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